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We present a mass determination for the transiting super-Earth Cancri e 
based on nearly 700 precise radial velocity (RV) measurements. This exten- 
sive RV data set consists of data collected by the McDonald Observatory planet 
search and published data from Lick and Keck observatories (Fischer et al. 2008). 
We obtained 212 RV measurements with the TuU Coude Spectrograph at the 
Harlan J. Smith 2.7m Telescope and combined them with a new Doppler re- 
duction of the 131 spectra that we have taken in 2003-2004 with the High- 
Resolution-Spectrograph (HRS) at the Hobby-Eberly Telescope (HET) for the 
original discovery of Cancri e. Using this large data set we obtain a 5-planet 
Keplerian orbital solution for the system and measure an RV semi-amplitude of 
X = 6.29 ± 0.21 ms-^ for Cnc e and determine a mass of 8.37±0.38 M®. The 
uncertainty in mass is thus less than 5%. This planet was previously found to 
transit its parent star (Winn et al. 2011, Demory et al. 2011), which allowed them 
to estimate its radius. Combined with the latest radius estimate from Gillon et 
al. (2012), we obtain a mean density of p = 4.50 ± 0.20 gcm^^. The location 
of p^ Cnc e in the mass-radius diagram suggests that the planet contains a sig- 
nificant amount of volitales, possibly a water-rich envelope surrounding a rocky 
core. 

Subject headings: planetary system — stars: individual (p^ Cancri, 55 Cancri, 
HR 3522, HD 75732) — techniques: radial velocities 



1. Introduction 

The p^ Cancri planetary system is one of the most interesting, nearby multi-planet 
systems that was discovered and extensively studied by the radial velocity (RV) technique. 
The parent star, p^ Cancri (= HR 3522, HD 75732, 55 Cancri), is a F = 5.95 G8V (Montes 
et al. 2001) star, located at a distance of 12.3 ± 0.1 pc, based on the Hipparcos parallax 
of 81.03 ± 0.75 mas (Van Leeuwen 2007). It is the primary of a wide visual binary with a 
projected separation of ^ 1065 AU (Mugrauer et al. 2006). The first planet in this system 
(p^ Cnc h, P — 14.65 d) was found by Butler et al. (1996), based on Lick Observatory RV 
data. Six years later, Marcy et al. (2002) presented evidence for two more giant planets 



^Based partly on observations obtained with the Hobby-Eberly Telescope, which is a joint project of the 
University of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximilians- 
Universitat Miinchen, and Georg-August-Universitat Gottingen. 



-3- 



orbiting this star: Cnc c with a period of 44 d and a very long period planet at an 
orbital separation of 5.5 AU, Cnc d. An intensive RV campaign that we carried out using 
the Hobby- Eberly Telescope (HET), revealed a short-periodic signal that was, at that time, 
thought to be one of the first discoveries of a hot Neptune with a minimum mass of ~ 17 
and an orbital period of 2.8 days (McArthur et al. 2004, hereafter Mc04), raising the total 
number of detected planets in this system to four. The planet count was further increased 
by Fischer et al. (2008, hereafter F08), who presented evidence for a fifth planet with an 
orbital period of 260 days. 

Dawson & Fabrycky (2010, hereafter DFIO) re-analyzed the published RV time-series 
data for this star and claimed that the 2.8 day period of Cnc e is an alias and that the true 
period of this companion is just 0.7365 days. This shorter period also led to a reduction of 
the minimum mass of p^ Cnc e to around 8 M^, moving the planet from the Neptune mass 
range into the super-Earth mass bin. DFIO pointed out the high a-priori transit probability 
of 25% for a planet and motivated highly precise photometric observations to search for the 
planetary transit signal. This transit search was indeed successful, using two different space 
telescopes, MOST and warm Spitzer. Winn et al. (2011, hereafter Wll), and Demory et 
al. (2011, hereafter Dll) practically simultaneously announced the detection of the transit 
signal of this super-Earth. The transit thus confirmed the shorter orbital period for p^ Cnc e, 
as suggested by DFIO. Dll measure a planetary radius of 2.1 ±0.17 i?© and a mean density 
of 4.8 ± 1.3gcm~^, while Wll reported a similar radius of 2.0 ± 0.14 i?© and shghtly higher 
mean density of 5.9ti;f gcm~^. This places p^ Cnc e in the mass-radius diagram between 
the region of high-density, rocky planets like CoRoT-7b (Hatzes et al. 2011) and Kepler- 10b 
(Batalha et al. 2011) and planets with a significant amount of volatiles, so-called "Mini- 
Neptunes", like GJ 1214b (Charbonneau et al. 2009) and the Kepler-11 planets (Lissauer et 
al. 2011). An improved radius determination of 2.17 ± 0.10 i?® was presented by Gillon et 
al. (2012, hereafter G12) by combining the Spitzer with the MOST photometry. Even more 
recently, Demory et al. (2012) used Warm Spitzer 4.5 pm observations of occultations of 
p^ Cnc e to detect its thermal emission. Clearly, this nearby transiting super-Earth planet 
around a bright star offers an abundant variety of very interesting follow-up observations 
that will allow a detailed characterization of this exoplanet. 

In our paper we will focus on the mass determination for p^ Cnc e based on hundreds 
of precise RV measurements. The paper is structured as follows: we first describe the 
observations of p^ Cnc at McDonald Observatory, the second section contains a description 
of the multi-planet orbital fit that we have performed, and the third section discusses our 
results, in particular the precise mass, for p^ Cnc e. 
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2. McDonald Observatory Observations 

We observed Cancri as part of the long-term Doppler exoplanet survey at the Harlan 
J. Smith 2.7m Telescope (HJST) (e.g. Cochran et al. 1997, Robertson et al. 2012) beginning 
in May 1999. For all observations we used the TuU Coude Spectrograph (TuU et al. 1995) 
in combination with a 1.2" slit that yields a resolving power oi R = 60,000 with 2-pixel 
sampling. The star light passes through a temperature controlled iodine vapor (I2) cell that 
is mounted in front of the entrance slit. The l2-cell provides a precise wavelength calibration 
and allows the reconstruction of the instrumental profile at the time of observation. Precise 
differential RVs are computed from each spectrum using our Austral Doppler code (Endl et 
al. 2000). 

As a consequence of the claim of DFIO that the true period of the inner planet is just 
0.7365 days, we changed our observing strategy for this star to test this prediction. While 
Cancri was previously observed with a cadence of once a month or less, we increased the 
cadence dramatically in 2010. To allow a good sampling of such a short period we typically 
observed p^ Cancri several times per month and often up to three times per night. Through 
June 2012 we collected a total of 212 single exposures, expanding the total time baseline of 
the HJST RV data to 4728 days or almost 13 years. 

We combine these RV measurements with an improved Doppler-reduction of the 131 
spectra obtained with the High-Resolution-Spectrograph (HRS) (TuU et al. 1995) at the 
Hobby-Eberly Telescope (HET) (Ramsey et al. 1998) over the time span of 190 days from 
October 2003 to April 2004. 

The 343 RV results from HJST/TuU and HET/HRS are displayed in Figured] and listed 
in Tabled] The large scatter in the RV results shown in Figured] is almost entirely due to 
Keplerian motion caused by the planetary system. 

3. Keplerian Orbital Solutions & Mass Determination 

We used GaussFit (Jefferys et al. 1988) to perform a simultaneous fit to our RV data 
of the known four giant planets in this system with periods of 14.7, 43, 260, and ~ 4700 
days. We do not add an additional noise term, often called "jitter", to the uncertainties of 
the RV measurements. Jitter is sometimes used to account for instrinsic stellar variability, 
algorithmical noise and/or residual instrumental systematic effects (Wright 2005). However, 
as the physical cause is poorly understood we decided not to add any jitter to the RV data. 
In principle, also the presence of additional planets can be a cause of larger residual scatter. 
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The 4-planet solution yields periods of Pl=14.6514 ± 0.0002 d, P2=43.0 ± 0.01 d, 
P3=261.8 ± 0.3 d and P4=4660 ± 37 d. The total rms-scatter of the 343 RV measurements 
around the 4-planet model fit is 7.5 ms~^ (HJST/Tull: 7.7ms~^, HET/HRS: 7.2ms~^). The 
reduced of this fit is 7.6. The biggest difference to F08 is the shorter period for the outer- 
most gas giant. However, the F08 period of 5218 ± 230 d is longer than the time baseline of 
our observations which probably explains why we do not find this longer period. We tested 
this by including the Lick data in the 4-planet model, which yield 4929 ±45 d for the best-fit 
period of planet d. 

We then searched the McDonald Observatory RV residuals from the 4-planet fit for a 
periodic signal due to the inner-most companion. The Lomb-Scargle periodogram of the 
residuals is displayed in Fig.[2]and has two prominent and highly significant peaks. The less 
significant peak with a Scargle-power of 38 is at a period of 2.8 days, the original period of 
planet e announced by M04. The shorter period of 0.7365 days is clearly the stronger peak 
with a power of 44.6. 

For all subsequent orbital fits we fixed the eccentricity of Cnc e to because nearly 
circular orbits are expected based on tidal damping (DFIO) and from classic secular theory 
(Van Laerhoven & Greenberg 2012). The observed occulation phase from D12 also sets an 
upper limit of e < 0.06 at the 3cr level, consistent with this expectation. 

In order to derive a precise companion mass from our measured K amplitude we also 
need a precise mass for the host star. We adopt the latest value of 0.905 ± 0.015 Mq from 
von Braun et al. (2011). 

We then fit a 5- planet model to our RV data with P = 0.7365 d for the inner- most 
companion. The 5-planet Keplerian fit yields a reduced of 3.8 and the residual rms-scatter 
around the fit is now reduced to 5.3ms^^ (HJST/Tull: 5.5ms~^, HET/HRS: 5.2ms~^). For 
Cnc e we obtain an RV semi-amplitude of K = 6.01 ± 0.38 ms^^. The errors produced 
by GaussFit are generated from a maximum likelihood estimation that is an approximation 
to a Bayesian maximum a posteriori estimator with a fiat prior (Jefferys 1990). Using the 
McDonald data alone derive a minimum mass of 7.95 ± 0.57 for p^ Cnc e. The stellar 
mass uncertainty of 1.7% and the K amplitude error of 6.2% hence lead to a planetary mass 
(l-cr) uncertainty of 7.2%. 

To further improve the precision of this mass determination we add the 70 published 
RV measurements from Keck Observatory from F08. The Keck data have a time baseline 
of 2037 days. We performed a new 5-planet fit including the abritrary RV zero-point of the 
Keck data set as a free parameter. The residual scatter of the 413 RV points is 5.5 ms~^ 
(HJST/Tull: 6.1ms-\ HET/HRS: 5.6ms-\ and Keck/HIRES: 3.1 ms^^). Again, we do 
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not detect any significant signal in the RV residuals. The HJST+HET+Keck 5-planet fit 
yields K = 6.22 ± 0.24 ms~^ and a minimum mass of 8.20 ± 0.41 M^. The addition of the 
Keck data increases the precision of the planet's mass determination to 5%. 

In a final step we also include the published long-term (6642 days) RV data (250 mea- 
surements) from Lick Observatory from F08. This raises the total number of precise RV 
measurements to 663 and the total time coverage to 8476 days or 23.2 years. The fi- 
nal 5-planet model using HJST+HET+Keck -|-Lick has a residual rms-scatter of 6.3 ms~^ 
(HJST/Tull: 6.1ms-S HET/HRS: 5.6ms-\ Keck/HIRES: 3.3ms-i and Lick/Hamilton: 
7.6 ms~^) . The Lomb-Scargle periodogram of the residuals shows no significant peak rising 
over the noise (see Figure[3]). Despite over 600 precise RV measurements we do not have any 
indication yet for a 6th planet in the system. 

Including the Lick data changes the value for the semi-amplitude of Cnc e slightly to 
K = 6.29 ± 0.21 ms~^ and the minimum mass to 8.30 ± 0.38 M®. By using the McDonald 
data together with the Keck and Lick RVs we thus arrive at a mass uncertainty of 4.6%. 
The phase-folded RV data are displayed in FigJH along with the best-fit Keplerian solution 
for Cnc e after subtracting the four other planets. Table[2] summarizes the parameters we 
obtained for p^ Cnc e from this full 5-planet solution. The orbital parameters for the four 
giant planets in the p^ Cnc system are listed in TableO Since we use the entire RV data, 
that were accumulated by the McDonald Observatory and California planet search, these 
orbital parameters should also be the best values to the 5-planet system. 



4. Discussion 

With the inclination i of 82.5 degrees measured by G12 from the photometric transit, 
we derive a planetary mass of 8.37 + 0.38 for p^ Cnc e. This value is slightly higher than 
reported by Dll (7.81 + 0.56 M®) and slightly lower than the value of Wll (8.63 + 0.35 M®). 
The Wll mass and its error were taken from the orbital solutions of DFIO. With our 
approach, we could not reproduce the small error in K of only 0.2 ms~^ presented by DFIO 
using the Keck and Lick RV data alone (we obtain an uncertainty in K of 0.33 ms^^). 

Our mass estimate is based on twice the amount of RV data and we also benefit from 
the improved stellar mass determination of von Braun et al. (2011). We therefore regard this 
new mass determination as the current best value for p^ Cnc e. The uncertainty in mass is 
driven by the uncertainty in K and in second order coupled to the uncertainty in the mass of 
the host star. We are thus only limited by the precision of the radial velocity measurements 
and the quality of the stellar parameters of p^ Cnc. 
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Figure|5] shows Cnc e in the mass-radius diagram compared to models for internal 
composition of small and low-mass planets from Seager et al. (2007) and three other transiting 
planets with well determined masses and radii, Cnc e has the smallest area of uncertainty 
based on its errors in mass and radius (Kepler- 10b has a smaller radius error but a larger 
uncertainty in mass). As noted by other authors (e.g. Gil), p^ Cnc e requires a significant 
amount of volatiles to explain its location in this diagram. It is located significantly above 
the model curves for purely rocky planets and approaches the zone of "mini-Neptunes" . Still, 
it does not require a large H/He envelope, but its mass and radius rather suggest a water-rich 
envelope around a rocky core. 

Kaib, Raymond & Duncan (2011) suggest that the p^ Cnc planetary system is coplanar 
but misaligned with its host star spin axis due to the perturbations of the secondary star. In 
principle, this can be tested (at least for p^ Cnc e) by observing and measuring the Rossiter- 
McLaughlin (RM) effect. However, the expected amplitude of the RM effect for p^ Cnc e is 
only 0.5 ms~^, and while RV measurements were obtained during transit by coincidence, a 
signal of this small amplitude is clearly undetectable by the current data. Future extreme 
precision RV measurements of several transits (by e.g. the upgraded HET/HRS or HARPS- 
North) might allow us to measure the spin-orbit misalignment for this planet. 

As mentioned before, despite the large quantity of precise RV measurements we did 
not detect any significant residual signal that could indicate a sixth planet in the system. 
The most interesting peak is near 131 days, as the window function is clean at this period 
value, and it would be at an orbital separation between the 43 d and the 261 d planet. But 
this peak has only a modest power and is not statistically signifcant. However, there is a 
rapidly increasing "treasure trove" of precise RV measurements for this system, with our 
paper adding over 300 RV points to the published sample. This should allow in the future to 
achieve sensitivity for more planets, either with lower mass or at longer orbital periods and 
especially in the habitable zone and in the currently large empty region between the inner 4 
planets and the distant outer planet at a 5 AU. The p^ Cnc multi-planetary system likely 
has more exciting discoveries waiting to be made. 

We thank our referee Kaspar von Braun for his thoughtful comments that helped to 
improve the manuscript. This material is based on work supported by the National Aero- 
nautics and Space Administration under Grant NNX09AB30G through the Origins of Solar 
Systems program. P. Robertson is supported by a University of Texas at Austin Con- 
tinuing Fellowship. The Hobby-Eberly Telescope (HET) is a joint project of the Univer- 
sity of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig- 
Maximilians-Universitat Miinchen, and Georg-August-Universitat Gottingen. The HET is 
named in honor of its principal benefactors, William P. Hobby and Robert E. Eberly. We 
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Tabic 1. Precise Radial Velocity Results 



BJD 



dRV(ms ^) err(ms ^) Telescope/Spectrograph 



51326.629453 
51503.889809 
51529.976108 



-22567.53 
-22568.26 
-22626.99 



11.8 
2.7 
1.3 



HJST/TuU 
HJST/TuU 
HJST/TuU 



Table 2. Keplerian Orbital Solution for the super-Earth planet Cnc e. 



parameter 


value 


notes 


Period (days) 


0.736546 ±0.000003 




K (ms-^) 


6.30 ±0.21 




e 


0.0 


(fixed) 


u (deg) 


90 


(fixed) 


Ttransit(<^ D) 


55568.011 ± 0.008 


55568.03469 


RMS (ms-i) 


6.09 


HJST/Tull (7 = -22574.1 ± 0.7 m s-i ) 


RMS (ms-i) 


5.62 


HET/HRS (7 = 28394.0 ± 0.08 ms^^ ) 


RMS (ms-i) 


3.25 


Kcck/HIRES (7 = 14.4 ± 0.5 ms-i ) 


RMS (ms-^) 


7.55 


Lick/Hamilton (7 = 14.2 ± 0.5 ms~^ ) 


mass (M®) 


8.37 ±0.38 


i = 82.5 


p (gcm^^) 


4.50 ±0.20 





Table 3. Keplerian Orbital Solutions for the 4 giant planets in the p^ Cnc system. 



parameter planet b planet c planet f planet d 



Period (days) 14.651 ± 0.0001 44.38 ± 0.007 261.2 ± 0.4 4909 ± 30 

K (ms-i) 71.11 ±0.24 10.12 ± 0.23 6.2 ± 0.3 45.2 ± 0.4 

e 0.004 ±0.003 0.07 ±0.02 0.32 ±0.05 0.02 ±0.008 

oj (deg) 110 ± 54 356 ± 22 139 ±8 254 ± 32 

Tpcr (BJD) 53035.0 ±2.2 53083 ± 3 51878 ± 5 53490 ± 437 

Msinz (Mjup) 0.80 ±0.012 0.165 ± 0.0054 0.172 ±0.008 3.53 ± 0.08 

a (AU) 0.1134 ±0.0006 0.237 ±0.0013 0.77 ±0.005 5.47 ± 0.06 
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Fig. 1. — 13 years of precise RV measurements of Cnc by the McDonald Observatory 
planet search. The filled (blue) circles are the HJST/Tull RV data and the filled (red) 
triangles represent the HET/HRS results. The two dense clusters of points are the high 
cadence campaigns we performed in 2003/2004 with the HET/HRS and starting in 2010 
with the HJST/Tull. We subtracted from both data sets the arbitrary, best-fit RV zero- 
points (HJST: 7 = — 22574. lms~^, HET: 7 = 28394 ms~^) from the joint Keplerian orbit 
fit (see section[3]). The large scatter of these data are almost entirely due to the Keplerian 
motion due to the planets in the Cnc system. 
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Fig. 2. — Lomb-Scargle periodogram of the 343 RV residuals after subtracting the orbits of 
the 4 known giant planets. Clearly, the shorter period at 0.74 d, as suggested by DFIO, is a 
stronger peak than the alias at 2.8 d. 
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Fig. 3. — Lomb-Scargle periodogram of the 663 RV residuals after subtracting the orbits 
of the 5 planets in the Cnc system. There is no significant signal apparent that could 
indicate the presence of a 6th planet in the system. 



- 14 - 




Fig. 4. — The 663 RV measurements of Cnc from HJST (filled circles), HET (diamonds), 
Keck (boxes) and Lick (triangles), after subtracting the four RV-orbits of the giant planets 
and phased to the transit ephemeris of Dll (transit occurs at phase 1, indicated by the 
vertical dashed line). Our best-fit Keplerian orbit for the super-Earth is shown as solid line. 
The width of this line represents the 1-a uncertainty of 0.21 ms~^ in the RV semi-amplitude 
K. The residual scatter around this orbit is 6.3 ms~^ (some of the worst outliers are outside 
the RV range shown here). 
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Fig. 5. — Cnc e in the mass-radius diagram and compared to other transiting super- 
Earths. The radius value was taken from G12 and the mass is from this work. The parameters 
for CoRoT-7b were adopted from Hatzes et al. (2011) and Kepler-lOb from Batalha et 
al.(2011). The radius value of GJ 1214b represent an average of the two estimates presented 
in Carter et al. (2011). The models for the internal composition of planets in this mass and 
radius range are taken from Seager et al. (2007). 



